qHTSWaterfall-vignette

library(qHTSWaterfall)

# Specify a comma separateed input file. In this example, one is included as 
# part of the qHTSPackage. If using your own file, provide the file path here.

# This example features active compounds from this study:
# Martinez NJ, Braisted JC, Dranchak PK, Moran JJ, Larson H, Queme B, Pak E, Dutra A, Rai G, 
# Cheng KC, Svaren J, Inglese J. Genome-Edited Coincidence and PMP22-HiBiT Fusion Reporter Cell 
# Lines Enable an Artifact-Suppressive Quantitative High-Throughput Screening Strategy for PMP22
# Gene-Dosage Disorder Drug Discovery. ACS Pharmacol Transl Sci. 2021 Jun 10;4(4):1422-1436.
# doi: 10.1021/acsptsci.1c00110. PMID: 34423274; PMCID: PMC8369676.


inputFile <- system.file("extdata", "NCATS_CMT1A_PMP22_Follow_Up.csv", package="qHTSWaterfall")

# Non-PubChemm input requires a list of log-molar concentrations, from low to high #concentration.
logConc <- c(
-9.011761791,
-8.534640544,
-8.057519287,
-7.580398031,
-7.103276777,
-6.626155522,
-6.149034267,
-5.671913012,
-5.194791758,
-4.717670503,
-4.240549248
)

# the main method call with parameters. Use 'help(plotWaterfall' for parameter descriptions.)
# Note: The displayed resulting plot shown in this vignette page is a low resolution static image.
# The method will return a high resolution 3D rendering (from the rgl package) that can be resized
# and rotated.
responseAxisConfig = list(min=-115, max=50, tickSizeVal=25, firstTick=-100)
concAxisConfig = list(min=-9.2, max=-3.9, tickSizeVal=1.0, firstTick=-4)
planeColors = list(basePlaneColor="#b8b6b6", rightPlaneColor="#999494", leftPlaneColor="#6e6868")

qHTSWaterfall::plotWaterfall(
  inputFile,
  activityReadouts = c("nluc","fluc"),
  pointColors = c("royalblue3", "darkgreen"),
  curveColors = c("royalblue3", "darkgreen"),
  inactiveColor = "gray",
  pointSize = 1.0,
  alpha = 1,
  plotInactivePoints = F,
  curveResolution = 250,
  plotAspectRatio = c(1.5,1,3),
  lineWeight = 1.5,
  fileFormat = 'ncats_qhts',
  planeColors = planeColors,
  showCurveNumberLabels = T,
  concAxisConfig = concAxisConfig,
  responseAxisConfig = responseAxisConfig
)
#> [1] "starting plotWaterfall method.."
#> Warning in qHTSWaterfall::plotWaterfall(inputFile, activityReadouts =
#> c("nluc", : NAs introduced by coercion

#> Warning in qHTSWaterfall::plotWaterfall(inputFile, activityReadouts =
#> c("nluc", : NAs introduced by coercion

#> Warning in qHTSWaterfall::plotWaterfall(inputFile, activityReadouts =
#> c("nluc", : NAs introduced by coercion

#> Warning in qHTSWaterfall::plotWaterfall(inputFile, activityReadouts =
#> c("nluc", : NAs introduced by coercion